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ABSTRACT 

We present a catalogue of X-ray selected galaxy clusters and groups as a first release of the 2XMMi/SDSS Galaxy Cluster Survey. The 
survey is a search for galaxy clusters detected serendipitously in observations with XMM-Newton in the footprint of the Sloan Digital 
Sky Survey (SDSS). The main aims of the survey are to identify new X-ray galaxy clusters, investigate their X-ray scaling relations, 
identify distant cluster candidates, and study the correlation of the X-ray and optical properties. In this paper, we describe the basic 
strategy to identify and characterize the X-ray cluster candidates that currently comprise 1 180 objects selected from the second XMM- 
Newton serendipitous source catalogue (2XMMi-DR3). Cross-correlation of the initial catalogue with recently published optically 
selected SDSS galaxy cluster catalogues yields photometric redshifts for 275 objects. Of these, 182 clusters have at least one member 
with a spectroscopic redshift from existing public data (SDSS-DR8). We developed an automated method to reprocess the XMM- 
Newton X-ray observations, determine the optimum source extraction radius, generate source and background spectra, and derive the 
temperatures and luminosities of the optically confirmed clusters. Here we present the X-ray properties of the first cluster sample, 
which comprises 175 clusters, among which 139 objects are new X-ray discoveries while the others were previously known as X- 
ray sources. For each cluster, the catalogue provides: two identifiers, coordinates, temperature, flux [0.5-2] keV, luminosity [0.5-2] 
keV extracted from an optimum aperture, bolometric luminosity L500, total mass M500, radius R^oo, and the optical properties of the 
counterpart. The first cluster sample from the survey covers a wide range of redshifts from 0.09 to 0.61, bolometric luminosities 
L500 = 1-9 X 10**^ - 1.2 X 10*' erg s"', and masses M500 - 2.3 x 10'"* - 4.9 x 10'* Mq. We extend the relation between the X-ray 
bolometric luminosity L500 and the X-ray temperature towards significantly lower T and L and still find that the slope of the linear 
L-T relation is consistent with values published for high luminosities. 
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1. Introduction 

Galaxy clusters are the most visible tracers of large-scale struc- 
ture. They occupy very massive dark matter halos and are obser- 
vationally accessible by a wide range of means. Their locations 
are found to corresponding to large numbers of tightly clustered 
galaxies, pools of hot X-ray emitting gas, and relatively strong 
features in the gravitational lensing shear field. Precise observa- 
tions of large numbers of clusters provide an important tool for 
testing our understanding of cosmology and structure formation. 
Clusters are also interesting laboratories for the study of galaxy 
evolution under the influence of extreme environments ( |Koester| 
let al.|2007| l. 

The baryonic matter of the clusters is found in two forms: 
first, individual galaxies within the cluster, which are most effec- 
tively studied through optical and NIR photometric and spectro- 
scopic surveys; and second, a hot, ionized intra-cluster medium 
(ICM), which can be studied by means of its X-ray emission and 
the Sunyaev-Zeldovich (SZ) effect ( jSunyaev & Z eldovi ch|1972] 
1980| l. The detection of clusters using SZ effect is a fairly new 
and highly promising technique for which tremendous progress 
has been made in finding high redshift clusters and measur 
ing the total cluster mass (e.g. Planck collaboration et al.|2011 



Vanderhnde et al.|2010[[Marriage et al.|2010[ ) 

The X-ray selection of clusters has several advantages for 
cosmological surveys: the observable X-ray luminosity and tem- 



perature of a cluster is tightly correlated with its total mass, 
which is indeed its most fundamental parameter (Reiprich & 
Bohringer 2002[ l. These relations provide the ability to mea- 
sure both the mass function (Boh ringer et al.| [2002) and power 
spectrum (Schuecker et al. 2003 ), which directly probe the cos- 
mological models. Since the cluster X-ray emission is strongly 
peaked on the dense cluster core. X-ray selection is less affected 
by projection effects than optical surveys and clusters can be 
identified efficiently over a wide redshift range. 

[Vlany clusters have been found in X-ray observations with 
Uhuru, HEAO-1, Ariel- V, Einstein, and EXOSAT, which have 
allowed a more accurate characterization of their physical pro- 
prieties (for a review see|Rosati et al.|(|2002| )). The ROSAT All 
Sky Survey (RASS, |Voges et al.||1999| l and the deep pointed 
observations have led to the discovery of hundreds of clus- 
ters. In ROSAT observations, 1743 clusters have been identi- 
fied, which are compiled in a meta-catalogue called MCXC 
by [Piffaretti et~aL| ( |2010[ ). The MCXC catalogue is based on 
pubhshed RASS-based (NORAS, REFLEEX, BCS, SGR NEP, 
MACS, and CIZA) and serendipitous (160D, 400D, SHARC, 
WARPS, and EMSS) cluster catalogues. 

The cuiTent generation of X-ray satellites XMM-Newton, 
Chandra, and Suzaku have provided follow-up observations of 
statistical samples of ROSAT clusters for cosmological studies 
( |Vikhlinin et al.|[2009| l and detailed information on the struc- 
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tural proprieties of the cluster population (e.g. |Vikhlinin et al.| 
2006[|Pratt et al.|2010[|Arnaud et al.|2010[ ). Several projects are 



ongoing to detect new clusters of galaxies from XMM-Newton 



and Chandra observations (e.g. the XSC (Ro mer et al.| 
XDCP dFassb ender et al.|2007|), XM M-LSS dPierre et al 



2001 



2006: 



COSM OS ([Finogueno vet al.||2007[), SXD S HFinoguenovlFa 
[20T0I 1, and ChaMP (Bai'kho use et al.|2 006)). 

In this paper, we present the 2XMMi/SDSS galaxy cluster 
survey, a search for galaxy clusters based on extended sources in 
the 2XMMi catalogue ( jWatson et al.|2009 ) in the field of view 
the Sloan Digital Sky Survey (SDSS). The main aim of the sur- 
vey is to build a large catalogue of new X-ray clusters in the sky 
coverage of SDSS. The catalogue will allow us to investigate the 
correlation between the X-ray and optical properties of the clus- 
ters. One of the long term goals of the project is to improve the 
X-ray scaling relations, and to prepare for the eROSITA cluster 
surveys, a mid term goal is the selection of the cluster candidates 
beyond the SDSS-limit for studies of the distant universe. Here 
we present a first cluster sample of the survey which comprises 
175 clusters found by cross-matching the 2XMMi sample with 
published SDSS based optical cluster catalogues. 

The paper is organized as follows: in Sect. 2 we describe 
the procedure of selecting the X-ray cluster candidates as well 
as their possible counterparts in SDSS data. In Sect. 3 we de- 
scribe the X-ray data the reduction and analysis of the optically 
confirmed clusters. The discussion of the results is described in 
Sect. 4. Section 5 concludes the paper The cosmological param- 
eters Qm = 0.3, Qa - 0.7 and //q = 70 km s"' Mpc"' were used 
throughout this paper. 



2. Sample construction 

We describe our basic strategy for identifying clusters among 
the extended X-ray sources in the 2XMMi catalogue. We then 
proceed by cross-matching the initial catalogue with those of 
optically selected galaxy clusters from the SDSS thus deriving 
a catalogue of X-ray selected and optically confirmed clusters 
with measured redshifts, whose X-ray properties are analysed in 
Sect. 3. 



2.1. X-ray cluster candidate list 

X-ray observations provide a robust method for the initial identi- 
fication of galaxy clusters as extended X-ray sources. A strategy 
to create a clean galaxy cluster sample is to construct a cata- 
logue of X-ray cluster candidates followed by optical observa- 
tions. XMM-Newton archival observations provide the basis for 
creating catalogues of serendipitously identified point-like and 
extended X-ray sources. The largest X-ray source catalogue ever 



produced is the second XMM-Newton source catalogue (Watson 
|et al.|2009| ). The latest edition of this catalogue is 2XMMi-DR3, 



which was released on 2010 April 28. The 2XMMi-DR3 cov- 
ers 504 deg^ and contains ~ 3 times as many discrete sources 
as either the ROSAT survey or pointed catalogues. The cata- 
logue contains 353191 X-ray source detections corresponding to 
262902 unique X-ray sources detected in 4953 XMM-Newton 
EPIC (European Photon Imaging Camera) observations made 
between 2000 February 3 and 2009 October 8. 

The 2XMMi-DR3 contains 30470 extended source detec- 
tions, which form the primary database for our study. This ini- 
tial sample contains a very significant number of spurious detec- 
tions caused by the clustering of unresolved point sources, edge 
effects, the shape of the PSF (point spread function) of the X- 



ray mirrors, large extended sources consisting of several minor 
sources, and other effects such as X-ray ghosts and similar 

We applied several selection steps to obtain a number of X- 
ray extended sources that were then visually inspected individ- 
ually. In our study, we considered only sources at high Galactic 
latitudes, \b\ > 20°, and discarded those that were flagged as 
spurious in the 2XMMi-DR3 catalogue by the screeners of the 
XMM-Newton SSC (Survey Science Centre). The source detec- 
tion pipeline used for the creation of the 2XMMi catalogues 
allows for a maximum core radius of extended sources of 80 
arcsec. Sources with extent parameters equalling that boundary 
were discarded, screening of a few examples shows that those 
sources are spurious or large extended sources (the targets of the 
observations) that were discarded anyhow. This initial selection 
reduced the number to 4027 detections. 

Since our main aim is the generation of a serendipitous clus- 
ter sample, we removed sources that were the targets of the 
XMM-Newton observation. We also discarded fields contain- 
ing large extended sources and selected only those fields within 
the footprint of SDSS, which left 1818 detections. After remov- 
ing multiple detections of the same extended sources (phenom- 
ena caused by either problematic source geometries in a sin- 
gle XMM-Newton observation or duplicate detections in a re- 
observed field), the catalogue was reduced to 1520 extended 
sources that were regarded as unique. 

This list still contains spurious detections for a number 
of reasons: (a) point-source confusion, (b) resolution of one 
asymmetric extended source into several symmetric extended 
sources, (c) the ill-known shape of the PSF leads to an excess 
of sources near bright point-sources, for both point-like and ex- 
tended sources, and (d) edge effects/low exposure times. To re- 
move the obvious spurious cases, we visually inspected the X- 
ray images of the initial 1520 detections using the FLIX upper 
limit serveiQ As a result, we were left with 1240 confirmed ex- 
tended X-ray sources. 

We then made use of the SDSS to remove additional non- 
cluster sources. We downloaded the XMM-Newton EPIC X-ray 
images from the XMM-Newton Science Archive (XSA: Arviset 
let al..2002 ) and created summed EPIC (PN+MOS1+MOS2) im- 
ages in the energy band 0.2 - 4.5 keV. Using these, we created 
smoothed X-ray contours, which were overlaid onto co-added r, 
i, and z-band SDSS images. Visual inspection of those optical 
multi-colour images with X-ray contours overlaid, allowed us to 
remove extended sources corresponding to nearby field galaxies, 
as well as those objects that are likely spurious detections. The 
resulting list which passes these selection criteria contains 1 180 
cluster candidates, about 75 percent of which are newly discov- 
ered. 

Figure [T] shows the X-ray-optical overlay of a new X-ray 
cluster, which has a counterpart in SDSS at photometric redshift 
- 0.4975 and a stellar mass centre indicated by the cross-hair as 
given by Szabo et al.|P011| (see Section 2.2). We use this clus- 
ter to illustrate the main steps of our analysis in the following 
sections. In Appendix A, Figures A.l to A.4 show the X-ray- 
optical overlays and the extracted X-ray spectra for four clusters 
illustrating results for various redshifts covered by our sample 
and different X-ray fluxes. 

About one quarter of the X-ray selected cluster candidates 
have no plausible optical counterpart. These are regarded as 
high-redshift candidates beyond the SDSS limit at z > 0.6, and 
suitable targets for dedicated optical/near-infrared follow-up ob- 
servations (see e.g. Lamer et al.|2008| l. 



http://www.ledas.ac.uk/flix/flix.html 
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Table 1. Main properties of the cluster catalogues with optically 
(SDSS-based) selected entries. The last two columns give the 
number of matching X-ray selected clusters individually and cu- 
mulatively. 



CLG 


Nr. 


Redshift 


SDSS 


X-ray 


Nr.CLG 


catalogue 


CLG 


range 




CLG(l') 


sample 


GMBCG 


55,000 


0.1-0.55 


DR7 


136 


123 


WHL 


39,688 


0.05 - 0.6 


DR6 


150 


72 


MaxBCG 


13,823 


0.1-0.3 


DR5 


54 


20 


AMF 


69,173 


0.045 - 0.78 


DR6 


127 


60 


Total 










275 




Fig. 1. The X-ray-optical overlay of the representative cluster 
2XMM J104421.8H-213029 at photometric redshift = 0.4975. 
The X-ray contours are overlaid on the SDSS co-added image 
obtained in r, /, and z-bands. The field of view is 4' x 4' centred 
on the X-ray cluster position. The cross-hair indicates the cluster 
stellar mass centre as given by |Szabo et al.| ( |20TT| ). 

2.2. The cross-matching with optical cluster catalogues 

The SDSS offers the opportunity to produce large galaxy- 
cluster catalogues. Several techniques were applied to identify 
likely clusters from multiband imaging and SDSS spectroscopy. 
We use those published catalogues to cross-identify common 
sources in our X-ray selected and those optical samples. All 
these optical catalogues give redshift information per cluster, 
which we use in the following to study the X-ray properties of 
our sources (zp indicates a photometric, Zs a spectroscopic red- 
shift). 

Table [T] lists the main properties of the optical cluster cata- 
logues that we used to confirm our X-ray selection. Below we 
provide a very brief description of each of these together with 
the acronym used by us: 

- GMBCG The Gaussian Mixture Brightest Cluster Galaxy 
catalogue ( |Hao et al.|2010) consists of more than 55,000 rich 
clusters across the redshift range 0.1 < Zp < 0.55 identified 
in SDSS-DR7. The galaxy clusters were detected by identi- 
fying the cluster red-sequence plus a brightest cluster galaxy 
(BCG). The cross-identification of X-ray cluster candidates 
with the GMBCG within a radius of 1 arcmin yields 136 
confirmed clusters. 

- WHL The catalogue of Wen, Han & Liu dWen et al.|2009] l 
consists of 39,668 clusters of galaxies drawn from SDSS- 
DR6 and covers the redshift range 0.05 < Zp < 0.6. A 
cluster was identified if more than eight member galaxies of 
Mr < 21 were found within a radius of 0.5 Mpc and within 
a photometric redshift interval Zp + 0.04(1 + Zp). We confirm 
150 X-ray clusters by cross-matching within 1 arcmin. 

- MaxBCG The max Brightest Cluster Galaxy catalogue 
( [Koester et al. 120071 lists 13,823 clusters in the redshift range 



0. 1 < < 0.3 from SDSS-DR5. The clusters were identified 
using maxBCG red-sequence technique, which uses the clus- 
tering of galaxies on the sky, in both magnitude and colour, 
to identify groups and clusters of bright E/SO red-sequence 
galaxies. The cross-match with our X-ray cluster candidate 
list reveals 54 clusters in common within a radius of one ar- 



AMF The Adaptive Matched Filter catalogue of Szabo et al. 
( 201 1 1 hsts 69, 173 likely galaxy clusters in the redshift range 



0.045 < Zp < 0.78 extracted from SDSS-DR6 using an adap- 
tive matched filter (AMF) cluster finder. The cross-match 
yields 127 confirmed X-ray galaxy clusters. 

In the AMF-catalogue, the cluster centre is given as the an- 
ticipated centre of the stellar mass of the cluster, while in the 
other three catalogues the cluster centre is the position of the 
brightest galaxy cluster (BCG). 

Many of our X-ray selected clusters have counterparts in sev- 
eral optical cluster catalogues within our chosen search radius 
of one arcmin. In these cases, we use the redshift of the optical 
counterpart, which has minimum spatial offset from the X-ray 
position. Table[T] lists in the second to last column the number 
of matching X-ray sources per optical catalog individually and 
in the last column the final number after removal of duplicate 
identifications. 

The unique optically confirmed X-ray cluster sample ob- 
tained by cross-matching with the four catalogues consists of 
275 objects having at least photometric redshifts. 

After cross-identification, we found 120 clusters of the opti- 
cally confirmed cluster sample with a spectroscopic redshift for 
the brightest galaxy cluster (BCG) from the published optical 
catalogues. Since the latest data release, SDSS DR8, provides 
more spectroscopic redshifts, we searched for additional spec- 
tra of BCGs and other member galaxies. We ran SDSS queries 
searching for galaxies with spectroscopic redshifts Zs(g) within 1 
Mpc from the X-ray centre. We considered a galaxy as a member 
of a cluster if \zp - Zs(g)\ < 0.05. 

The spectroscopic redshift of the cluster was calculated as 
the average redshift for the cluster galaxies with spectroscopic 
redshifts. The confirmed cluster sample with spectroscopic red- 
shifts for at least one galaxy includes 182 objects. Therefore, the 
unique optically confirmed X-ray cluster sample has the photo- 
metric redshifts for all of them, 120 spectroscopic redshifts for 
120 BCGs from the optical cluster catalogues, and 182 clusters 
with one or more members with spectroscopic redshifts from 
the SDSS database. Figure |2] shows the distribution of the clus- 
ter photometric redshift Zp, the distribution of spectroscopic red- 
shifts Zs of the BCGs as given in the various optical cluster cat- 
alogues, and the average spectroscopic redshift of the cluster 
members (which we refer to as the cluster spectroscopic redshift) 
for the confirmed cluster sample. Figure [3] shows the distribu- 
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Fig. 2. The distribution of the optical redshifts for the confirmed 
clusters sample. The distribution includes the cluster photomet- 
ric redshifts Zp (solid line) with a median 0.36, spectroscopic 
redshifts of the BCGs Zs (dashed line) with a median 0.3 from 
the optical cluster catalogues and the cluster spectroscopic red- 
shifts Zs (dotted line) with a median 0.3 from the SDSS data. 




Fig. 4. The relation between the photometric and spectroscopic 
redshifts of the confirmed cluster sample. 




100 150 200 250 300 
Optical/X-ray offset (kpc) 



10 20 30 40 50 60 

Number of cluster galaxies with spectroscopic redshifts 



Fig. 5. The distribution of the linear separation between the po- 
sition of the BCG and the X-ray cluster position. 



Fig. 3. The distribution of the spectroscopically confirmed clus- 
ter members per cluster. 



tion of the number of cluster galaxies that have a spectroscopic 
redshift in the SDSS database. The relation between the photo- 
metric and spectroscopic redshifts of the cluster sample is shown 
in Figure |4] Since this relation was found to be tight (where the 
Gaussian distribution of (zs - Zp) has cr = 0.02), we were able to 
rely on the photometric redshifts for the cluster with no spectro- 
scopic information. 

We used an angular separation of one arcmin to cross-match 
the X-ray cluster candidates with the optical cluster catalogues. 
The corresponding linear separation was calculated using the 
spectroscopic redshift, if available, or the photometric redshift. 
Figure |5] shows the distribution of the linear separation between 
the X-ray centre and the BCG position. For AMF clusters (60 ob- 
jects), we identified the BCGs of 40 systems within one arcmin 
and computed their off'sets, which are included in this afore- 
mentioned distribution. The BCGs were selected as the brightest 
galaxies with \zp - Zp(BCG cand.)\ < 0.05 among the three BCG 



candidates given for each AMF cluster published by Szabo et al. 
( 201 l| l. The other 20 AMF clusters are not included inFig.|5J be- 



cause their BCG is outside one arcmin. It is not always the case 
that the BCG lies exactly on the X-ray peak. lRykoff' et al.| ( [2008| 
model the optical/X-ray offset distribution by matching a sample 
of maxBCG clusters to known X-ray sources from the ROSAT 
survey. They found a large excess of X-ray clusters associated 
with the optical cluster centre. There is a tight core in which the 
BCG is within ~ 150 h"' kpc of any X-ray source, as well as a 
long tail extending to ~ 1500 h"' kpc. It is shown in Figure [s] 
that the majority of the confirmed sample have the BCG within 
a radius (~ 150 kpc), as well as a tail extending to 352 kpc that 



is consistent with the optical/X-ray offset distribution of Rykoff^ 
et al. (2008). 

We searched the Astronomical Database SIMBAD and the 
NASA/IPAC Extragalactic Database (NED) to check whether 
they had been identified and catalogued previously. We used a 
search radius of one arcmin. About 85 percent of the confirmed 
sample are new X-ray clusters, while the remainder had been 
previously studied using ROSAT, Chandra, or XMM-Newton 
data. 
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Fig. 6. The distribution of EPIC counts in [0.2-12] keV as given 
in the 2XMMi-DR3 catalogue for the confirmed cluster sample. 



3. X-ray data analysis 

The optically confirmed clusters have a wide range of source 
counts (EPIC counts in the broad band energy 0.2-12 keV from 
the 2XMMi-DR3 catalogue) from 66 to 28000 counts as shown 
in Figure 6] To analyse the X-ray data, we have to determine 
the optica redshifts, except for some candidates with more than 
1000 net photons for which it is possible to estimate the X-ray 
redshift (e.g. [ Lamer et alT][2008l |Yu et al.|[20n] l. In this paper, 
we use the cluster spectroscopic redshifts where available or the 
photometric redshifts that we obtained from the cross-matching 
as described in the previous section. 

The data reduction and analysis of the optically con- 
firmed sample was carried out using the XMM-Newton Science 
Analysis Software (SAS) version 10.0.0. 



3.1. Standard pipelines 

The raw XMM-Newton data were downloaded using the Archive 
Interoperability System (AIO), which provides access to the 
XMM-Newton Science Archive (XSA). The raw data were pro- 
vided in the form of a bundle of files known as observation data 
files (ODF), which contain uncalibrated event files, satellite atti- 
tude files, and calibration information. The main steps in the data 
reduction were: (i) the generation of calibrated event Usts for the 
EPIC (MOSl, M0S2, and PN) cameras using the latest calibra- 
tion data. This was done using the SAS packages cifbuild, 
odf ingest, epchain, and emchain. (ii) The creation of back- 
ground light curves to identify time intervals with poor quality 
data, (iii) The filtering of the EPIC event lists to exclude periods 
of high background flaring and bad events, (iv) To create a sky 
image of the filtered data set. The last three steps were performed 
using SAS packages evselect, tabgtigen, and xinmselect. 

3.2. Analysis of tlie sample 

We now describe the procedure to determine the source and 
background regions for each cluster, extract the source and back- 
ground spectra, fit the X-ray spectra, and finally measure the X- 
ray parameters (e.g. temperature, flux, and luminosity). As in- 
put to the task generating the X-ray spectra, we used the filtered 
event lists as described in the previous section. 




20 40 60 80 100 120 140 160 IBO 

Radius (arcsec) 



Fig.7. The radial profile of 2XMM J104421.8-H213029 
MOSl (green), MOS2(red), and PN (blue) images. The horizon- 
tal lines indicate the background values for MOSl, MOS2 and 
PN with the same colour as the profile. 



3.2.1. Optimum source extraction radius 

The most critical step in generating the cluster X-ray spectra is to 
determine the source extraction radius. We developed a method 
to optimize the signal-to-noise ratio (SNR) of the spectrum for 
each cluster. To calculate the extraction radius with the highest 
integrated SNR we created radial profiles of each cluster in the 
energy band 0.5 - 2.0 keV. Background sources, taken from the 
EPIC PPS source lists, were excluded and the profiles were ex- 
posure corrected using the EPIC exposure maps. Since we did 
not perform a new source detection run, the SNR was calculated 
as a function of radius taking into account the background levels 
as given in the 2XMMi catalogue. 

The radial profiles of the X-ray surface brightness of the rep- 
resentative cluster in MOSl, MOS2 and PN data are shown in 
Fig. [7] The background values of the cluster in the EPIC images 
are indicated by the horizontal line with the same colours as the 
profiles. Figure [8] shows the SNR profiles of the representative 
cluster in MOSl, MOS2, PN and EPIC (MOS1H-MOS2+PN) 
data as a function of the radius from the cluster centre. The op- 
timum extraction radius (72") is determined from the maximum 
value in the EPIC SNR plot, which is indicated by a point in 
Fig. [8] 

3.2.2. Spectral extraction 

The EPIC filtered event lists were used to extract the X-ray 
spectra of the cross-correlated X-ray optical cluster sample. The 
spectra of each cluster candidate were extracted from a region 
with an optimum extraction radius as described in the previous 
section. The background spectra were extracted from a circular 
annulus around the cluster with inner and outer radii equalling 
two and three times the optimum radius, respectively. Other 
unrelated nearby sources were masked and excluded from the 
source and background regions that were finally used to extract 
the X-ray spectra. Figure |9] shows the cluster and background 
regions, as well as the excluded regions of field sources for the 
representative cluster. The SAS task especget was used to gen- 
erate the cluster and background spectra and to create the re- 
sponse matrix files (redistribution matrix file (RMF) and ancil- 
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Fig. 8. The signal-to-noise ratio (SNR) profiles of 2XMM 
J104421.8+213029 in MOSl (green), MOS2 (red), PN (blue) 
and EPIC (MOS1+MOS2+PN) (cyan) data. The cluster opti- 
mum extraction radius (72") is corresponding to the highest 
SNR as indicated by a point in the EPIC SNR profile. 
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Fig. 10. The EPIC PN (black), MOSl (green) and MOS2 (red) 
data with the best-fit MEKAL model for the representative clus- 
ter. 




Fig. 9. The representative cluster extraction region is the inner 
circle with colour cyan. The background region is the annulus 
with white colour The excluding field sources are indicated by 
green circles. The field of view is 8' x 8' centred at the cluster 
position. 



lary response file (ARE)) required to perform the X-ray spectral 
fitting with XSPEC. 



3.2.3. Spectral fitting 

The photon counts of each cluster spectrum were grouped into 
bins with at least one count per bin before a fit of a spectral 
model was applied to the data using the Etools task grppha. The 
spectral fitting was carried out using XSPEC software version 



12.5.1 ( Arnaud|1996 1. Before executing the algorithm to fit the 
spectra, the Galactic HI column (nH) was derived from the HI 
map from the Leiden/Argentine/Bonn (LAB) survey (jKalberla 
et al. 2005). This parameter was fixed while fitting the X-ray 
spectrum. The redshift of the spectral model was fixed to the 
optical cluster redshift either the spectroscopic redshift for 182 
clusters or the photometric redshift for the remainder cluster 
sample. 

For each cluster, the available EPIC spectra were fitted si- 
multaneously. The employed fitting model was a multiplica- 
tion of a TBABS absorption model ( Wilms et al., 2000) and 
a single-temperature optically thin thermal plasma component 
(the MEKAL code in XSPEC terminology, I Mewe et aL]|1986| ) to 
model the X-ray plasma emission from the ICM. The metallicity 
was fixed at 0.4 Zo. This value is the mean of the metallici- 
ties of 95 galaxy clusters in the redshift range from 0. 1 to 0.6 
(the same redshift range of the confirmed sample) observed by 
Chandra ( Maughan et al.|2()08 l. The free parameters are the X- 
ray temperature and the spectral normalization. The fitting was 
done using the Cash statistic with one count per bin following the 



recommendation of Krumpe et al. ( 2008 1 for small count statis- 
tics. 

To avoid the fitting algorithm converging to a local mini- 
mum of the fitting statistics, we ran series of fits stepping from 
0.1 to 15 keV with a step size = 0.05 using the steppar com- 
mand within XSPEC. The cluster temperature, its flux in the 
[0.5-2] keV band, its X-ray luminosity in the [0.5-2] keV band, 
the bolometric luminosity, and the corresponding errors were de- 
rived from the best-fitting model. We assumed that the fractional 
error in the bolometric luminosity was the same as the fractional 
error in the aperture luminosity [0.5-2] keV (within an aperture 
defined by the optimum extraction radius). Figure 10 shows the 



fits to the EPIC (MOSl, MOS2, and PN) specti-a and the models 
for the representative cluster Figures A. 1 to A. 4 in the Appendix 
A show the fitted spectra of four clusters with different X-ray 
surface brightnesses and data qualities at different redshifts cov- 
ering the whole redshift range of the confirmed clusters. 
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Fig. 11. The distribution of the fractional errors in the X-ray tem- 
peratures (solid) and luminosities (dashed) derived from spectra 
extracted within the optimum aperture for the energy range 0.5-2 
keV of the first cluster sample. 



4. Analysis of a cluster sample with reliable X-ray 
parameters 

We analysed the X-ray data of the optically and X-ray con- 
firmed clusters to measure the average temperature of the hot 
ICM. We developed an optimal extraction method for the X- 
ray spectra maximising the SNR. The cluster spectra were fitted 
with absorbed thin thermal plasma emission models with pre- 
determined redshift and interstellar column density to determine 
the aperture X-ray temperature (Tap), flux (Fap) [0.5-2] keV, lu- 
minosity (L„p) [0.5-2] keV, and their errors. We accepted the 
measurements of Tap and L„p if the fractional errors were smaller 
than 0.5. About 80 percent of the confirmed clusters passed this 
fractional error filter. For these clusters, another visual screening 



of the spectral fits (Figure 10 1 and the X-ray images (Figure^ 



was done. When the spectral extraction of a given cluster was 
strongly affected by the exclusion of field sources within the ex- 
traction radius or a poor determination of the background spec- 
trum, it was also excluded from the final sample, which com- 
prises 175 clusters. For a fraction of 80 percent, this is the first 
X-ray detection and the first temperature measurement. 

Our subsequent presentation of our analysis and discussion 
refers to those 175 objects with reliable X-ray parameters. The 
distribution of the Tap and Lap [0.5-2] keV fractional errors for 
the first cluster sample is shown in Fig. [TT] It is clearly evident 
that the cluster luminosity is more tightly constrained than the 
temperature. For about 86 percent of the sample, the fractional 
errors are smaller than 0.25. Therefore, we estimated several 
physical parameters for each cluster based on the bolometric lu- 
minosity Lhoi within the optimal aperture. The median correction 
factor between aperture bolometric luminosities and aperture lu- 
minosities in the energy band [0.5-2] keV ( L^o/ / L„p ) was found 
to be 1.7. We assumed that the fractional error in Lhoi was iden- 
tical to that of Lap [0.5-2] keV. The estimated parameters are 
.^500 the radius at which the mean mass density is 500 times the 
critical density of the Universe (see Eq. 2) at the cluster red- 
shift, Lgoo the bolometric luminosity within Rsoo, and M500 the 
cluster mass within R^oo- We used an iterative procedure to es- 
timate the physical parameters using published Lsoo - T^qq and 
L500 - M500 relations ( [Pratt et al.||2009| their orthogonal fit for 



M500 with Mal mquist bias correction) . Ou r procedure is simila r 
to that used by 



Piff"aretti et al. 



(2010: and 



Suhada et al. 



which consists of the following steps: 
(i) We estimate M500 using the L - M relation 



(2010 1, 



500 



2x 10' 



1/2.08 



38 X 1044 erg s- 



(1) 



day value, h{z) - [nM(l +1^ + 



where h{z) is the Hubble constant normaUsed to its present- 

1/2 

. We approximate L500 

as the aperture bolometric luminosity Lhoi, which we correct 
in an iterative way. 
(ii) We compute Rsm 



Rm) = ^/3M5oo/4;r500p,(^), 



(2) 



where the critical density is pdz) = h(z)^3H^ l%nG . 
(iii) We compute the cluster temperature within R^qq using the 
L-T relation 



= 5keV (- 



h{z) ' Lt„i 



,1/3.35 



13 X 1044 erg i- 



(3) 



(iv) We calculate the core radius rcore and y6 using scaling rela- 
tions from Finogu enov et al.| ( |2007 1 



r^ore = 0.07 X 7^500 X 



xO.63 



keV 



) 



I ^ \l/3 



(4) 



(5) 



(v) We calculate the enclosed flux within /?5oo and the optimum 
aperture by extrapolating the /?-model. The ratio of the two 
fluxes is calculated, i.e. y = FjQo/F[,oi. 

(vi) We finally compute a corrected value of L500 - yx iboi- 

We then considered L500 as input for another iteration and all 
computed parameters were updated. We repeated this iterative 
procedure until converging to a final solution. At this stage, the 
L500, Msoo, and R^qo were determined. The median correction 
factor between extrapolated luminosities and aperture bolomet- 
ric luminosities (Lsoo/iboi) was 1.5. To calculate the errors in 
Eqs. 1 and 3, we included the measurement errors in the aper- 
ture bolometric luminosity Li,oi , the intrinsic scatter in the L-T 
and L- M relations, and the propagated errors caused by the un- 
certainty in their slopes and intercepts. For Eq. 4 and Eq. 5, we 
included only the propagated errors of their independent param- 
eters since their intrinsic scatter had not been published. Finally, 
all the measured errors were taken into account when comput- 
ing the errors in L500 and Mgoo in the last iteration. The errors in 
L500 and M500 were still underestimated because of the possible 
scatter in the relations in Eq. 4 and Eq. 5. 

We investigated the L-T relation in the first cluster sam- 
ple using Tap and Lsoo- Figure 12 shows the relation between 
L500 (corrected for the redshift evolution) and Tap (uncorrected 
for cooling flows). Here we assumed that Tap did not differ sig- 
nificantly from Tsoo and its error was derived from the spectral 
fits. The best-fit linear relation (solid line) derived from an or- 
thogonal distance regression fit (ODR) ( Boggs & Rogers] 1990| 
which takes into account measurement errors in both variables) 
between their logarithm, is 

log (hiz)-' Lsoo) = (0.57 ± 0.05) + (3.41 ± 0.15) log (Tap). (6) 
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Fig. 12. The relation between the X-ray bolometric luminosities 
L500 and aperture temperatures T^p of the first cluster sample. 
The solid line indicates the best fit of the sample using orthog- 
onal distance regression (ODR). The dashed line is the extrap- 
olated relation for REXCESS sample ( [Pratt et al.|2009| l using a 
BCES orthogonal fit. 



The best-fit power law relation derived from a BCES orthog- 
onal fit to the Lsoo - Tsoo relation published by Pratt et i 



for the REXCESS sample is plotted as the dashed hne in 
The ODR slope (present work), 3.41 + 0.15, is consistent with 



aL]( p009ll 
inFiiqi| 



the BCES orthogonal slope ( [Pratt et al.|2"009j of the REXCESS 
sample, 3.35 ± 0.32. In addition, the present slope is consistent 
with the BCES orthogonal slope (3.63 ± 0.27) of the L - T re- 
lation derived from a sample of 1 14 clusters (without excluding 
the core regions) observed with Chandra across a wide range of 
temperature (2 < kT < 16 keV) and redshift (0.1 < z < 1.3) by 
Maughan etaLl ( [201 l| l. 



We tested the corresponding uncertainty in the error budget 
of Lgoo caused by the above-mentioned unknown scatter in Eq. 
4 and 5: for example, a cr = +0.1 scatter in the /3 value results 
in a fractional error in L500 of 17%. If we take into account the 
newly estimated errors in L500 when fitting the L - T relation, 
the revised slope of 3.32 is within the error in the original slope 
as in Eq. 6 and still consistent with the published ones. 

Our sample represents cluster temperatures ranging from 
0.45 to 5.92 keV and values of bolometric luminosity in the L500 
range 1.9 x 10"*^ - 1.2 x 10"^^ erg s ' in a wide redshift range 
0.1 - 0.6. Most of the published L - T relations were derived 
from local cluster samples with temperatures higher than 2 keV. 
The current relation is derived for our sample, which includes 
clusters and groups with low temperatures and luminosities in a 
wide redshift range up to z = 0.6. The distribution of luminosity 
as a function of redshift is shown in Fig. 13 

Table 2, available at the CDS, represents the first cluster 
sample containing as many as 175 X-ray clusters. In addi- 
tion, the first cluster sample with the X-ray-optical overlay 
and fitted spectra for each cluster is publicly available from 
http://www.aip.de/groups/xray/XMM_SDSS_CLUSTERS In 
the catalogue, we provide the cluster identification number 
(detection Id, detid) and its name (lAUNAME) in (cols. [1] 
and [2]), the right ascension and declination of X-ray emission 
in equinox J2000.0 (cols. [3] and [4]), the XMM-Newton 
observation Id (obsid) (col. [5]), the optical redshift (col. [6]), 
the scale at the cluster redshift in kpc/" (col. [7]), the aperture 



Fig. 13. The distribution of the X-ray bolometric luminosities 
L500 as a function of optical redshifts of the first cluster sample. 



and /?5()o radii in kpc (col. [8] and [9]), the cluster aperture X-ray 
temperature Tap and its positive and negative errors in keV 
(cols. [10], [11] and [12], respectively), the aperture X-ray flux 
Fap [0.5-2] keV and its positive and negative errors in units 
of 10"'"* erg cm"^ s"' (cols. [13], [14] and [15], respectively), 
the aperture X-ray luminosity Lap [0.5-2] keV and its positive 
and negative errors in units of lO'*^ erg s"' (cols. [16], [17] 
and [18], respectively), the cluster bolometric luminosity L500 
and its error in units of 10"*^ erg s"' (cols. [19] and [20]), the 
cluster mass M^nn and its error in units of lO'"* Mo (cols. [21] 



and [22]), the Galactic HI column in units 10 cm" (col. [23]), 
the identification number of the cluster in optical catalogue 
(col. [24]), the BCG right ascension and declination in equinox 
J2000.0 (cols. [25] and [26]) although for AMF catalogue they 
represent the cluster stellar mass centre, the cluster photometric 
redshift (col. [27]), the average spectroscopic redshift of the 
cluster galaxies with available spectroscopic redshifts and 
their number (cols. [28] and [29]), the linear offset between the 
cluster X-ray position and the cluster optical position (col. [30]), 
the optical cluster catalogue names that identify the cluster 
(col. [31]) (Note: the optical parameters are extracted from 
the first one), and the alternative name of the X-ray clusters 
previously identified using ROSAT, Chandra, or XMM-Newton 
data and its reference in NED and SIMBAD databases (col. [32] 
and [33]). 



5. Summary and outlook 

We have presented the first sample of X-ray galaxy clusters 
from the 2XMMi-Newton/SDSS Galaxy Cluster Survey. The 
survey comprises 1180 cluster candidates selected as X-ray 
serendipitous sources from the second XMM-Newton serendipi- 
tous source catalogue (2XMMi-DR3) that had been observed by 
the SDSS. A quarter of the candidates are identified as distant 
cluster candidates beyond z - 0.6, because there is no appar- 
ent overdensity of galaxies in the corresponding SDSS images. 
Another quarter of the candidates had been previously identi- 
fied in optical cluster catalogues extracted from SDSS data. Our 
cross-correlation of the X-ray cluster candidates with four opti- 
cal cluster catalogues within a matching radius of one arcmin 
confirmed 275 clusters and provided us with the photometric 
redshifts for all of them and the spectroscopic redshifts for 120 
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BCGs. We extracted all available spectroscopic redshifts for the 
cluster members from recent SDSS data. Among the confirmed 
cluster sample, 182 clusters have spectroscopic redshifts for at 
least one galaxy member More than 80 percent of the con- 
firmed sample are newly identified X-ray clusters and the oth- 
ers had been previously identified using ROSAT, Chandra, or 
XMM-Newton data. We reduced and analysed the X-ray data 
of the confirmed sample in an automated way. The X-ray tem- 
perature, flux and luminosity of the confirmed sample and their 
errors were derived from spectral fitting. The analysed sample 
in the present work contains 175 X-ray galaxy clusters with ac- 
ceptable measurements of X-ray parameters (ie. with fractional 
errors smaller than 0.5) from reasonable quality fitting (139 ob- 
jects being newly discovered in X-rays). In addition, we derived 
the physical properties (Rsoq, L500 and M500) of the study sam- 
ple from an iterative procedure using the published scaling re- 
lations. The relation between the X-ray bolometric luminosity 
L500 and aperture temperature of the sample is investigated. The 
slope of the relation agrees with the slope of the same relation 
in the REXCESS sample ( Pratt et al.||2509| . The present rela- 
tion is derived from a large sample with low luminosities and 
temperatures across a wide redshift range 0.09 - 0.61. 

As one extension to this project, we intend to obtain SDSS 
photometric redshifts of all 2XMMi-DR3 X-ray cluster candi- 
dates that have been detected in SDSS imaging. This will signif- 
icantly increase the sample size and the identified fraction of the 
2XMMi cluster sample. Further improvements in the accuracy 
of the X-ray parameters for about 10 percent of the confirmed 
sample will be made by analysing repeated observations of those 
clusters. 
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Appendix A: Gallery 

We present a gallery of four galaxy clusters from the first clus- 
ter sample with different X-ray fluxes and data quality at differ- 
ent redshifts covering the whole redshift range of the sample. 
For each cluster. X-ray flux contours (0.2-4.5 keV) are overlaid 
on combined image from r, and z- SDSS images. The up- 
per panel in each figure shows the X-ray-optical overlays. The 
field of view is 4' x 4' centred on the X-ray cluster position. In 
each overlay, the cross-hair indicates the position of the bright- 
est cluster galaxy (BCG), while in Figure A.4 the cross-hair in- 
dicates the cluster stellar mass centre although but it is obvious 
that the BCG is located at the X-ray emission peak. In each fig- 
ure, the bottom panel shows the X-ray spectra (EPIC PN (black), 
MOSl (green), MOS2 (red)) and the best fitting MEKAL 
model. The full gallery of the first cluster sample is available 
at |http://www.aip.de/groups/xray/XMM_SDSS-CLUSTERS. 
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Fig.A.l. detid = 275409: 2XMMI J143929.0-H024605 at z, 
0.1447 (Fap [0.5 - 2] keV = 0.63 x lO"'"* erg cm^^ s"'). 
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as 1 2 5 

Energy (keV) 



Fig.A.4. detid = 097911: 2XMM J092545.5+305858 at Zp = 
0.5865 (Fap [0.5 - 2] keV = 7.59 x lO-'"* erg cm-^ g-'). 
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